Large Spin Torque in Topological Insulator/Ferromagnetic Metal Bilayers 
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While there are several proposals on spintronic applications of the surface state's helical spin 
texture in three-dimensional (3D) topological insulators (TI), none has been realized yet. Motivated 
by recent measurements of spin-transfer torque in heavy metal (HM)/ ferromagnetic-metal (FM) 
structures [1-3], we propose measuring spin-transfer torque in a TI/ FM bilayer with in-plane 
magnetic moment and in-plane current. Based on ARPES experiments on TI/FM bilayers and 
transport measurements on thin TI films, we focus on the role of two types of surface states: the 
topological Dirac surface state and an additional two-dimensional electron gas (2DEG) with Rashba 
spin-orbit coupling. We find that each surface state contributes to both an out-of-plane (field-like) 
and an in-plane torque. However, the torques due to the topological Dirac surface state and the 
Rashba-split 2DEG turn out to be in opposite directions. This will allow for the experimental 
assessment of the dominant source of spin-torque. Finally, we estimate the spin-torque efficiency in 
TI/FM bilayers to be comparable to that of the HM/FM system opening up exciting prospect for 
applications. 



Introduction - Theoretical proposals for harnessing the 
helical spin texture of topological surface states of strong 
3D TIs [4] so far have focused on Tl/ferromagnetic in- 
sulator (FI) hybrid structures [5-7]. In these (ideal) 
TI/FI bilayers, the topological surface state is well de- 
fined at the interface between the two insulators and 
carries all current running through the system. How- 
ever, room-temperature FIs are scarce. Though YIG 
is a rare room-temperature FI, it is not suitable for 
nanoscale fabrication. Furthermore, for certain applica- 
tions, such as a three-terminal device (magnetic mem- 
ory), a metallic magnetic material allows for a straight- 
forward (resistive) readout [2, 3]. Finally, the applica- 
bility of the ferromagnetic-resonance technique for mea- 
suring the spin-torque efficiency [1-3] makes a TI/FM 
bilayer an attractive system. Hence, we propose the ex- 
perimental investigation of spin torque due to a current 
running through a TI/FM bilayer [see Fig. 1(a)] and an- 
alyze here the role of the topological surface state and its 
helical spin texture. 

Current-driven manipulation of nanoscale magnetic 
devices holds great promise for next generation mem- 
ory devices. One such mechanism relies on a current 
through a (non-magnetic) material leading to a torque 
on an adjacent ferromagnetic metal through the spin-Hall 
effect. The figure of merit for such a device is the spin- 
Hall angle: a dimensionless ratio of spin-current density 
perpendicular to charge-current density. This quantity 
measures the efficiency of spin-torque generation driven 
by charge current. Such spin-Hall-cffcct based magne- 
tization switching has recently been demonstrated for 
Ta- and W-based devices with switching currents I < 1 
mA [2, 3]. However, for practical applications, even 
higher spin-torque efficiencies (lower switching currents) 
would be desirable. In this work, we propose using TIs for 




FIG. 1. (a) TI / FM bilayer we consider in this work with a 
current running in x direction and an in-plane magnetization 
M — Mm along the in-plane diagonal, (b) A schematic band- 
structure with a Dirac surface state and a Rashba-split 2DEG 
as seen in Refs.[9-12]. The corresponding spin structure at 
£f is shown in (c). 



the torque generation, utilizing the surface state's spin- 
momentum locking and predict a substantial spin-torque 
efficiency. 

We need to address two issues in order to establish a 
simplest theoretical model that captures salient realistic 
aspects of TI/FM bilayers: (1) the influence of the FM on 
the TI surface band structure, and (2) the possibility of 
TI bulk conduction. The first issue arises as the topolog- 
ical surface state is not necessarily protected when the TI 
is in contact with a metal [8]. However, ARPES results 
on TIs with thin FM layers on top show a bandstructure 
as sketched in Fig. 1(b): the topological surface state 
survives and additional (Rashba-split) surface bands are 
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formed due to band bending [9-12]. The second issue 
has to do with the chemical potential of most topologi- 
cal insulators, including Bi 2 Se3 the to date best studied 
material [13-16], lying in the conduction band. However, 
recent experiments showed that surface states dominate 
the transport for thin Bi2Se3 films of high-quality [17]. 
Based on these observations we will focus on the topo- 
logical (Dirac) state with its linear dispersion [thick solid 
line in Fig. 1(b) and (c)] and the Rashba-split 2DEG 
(dashed lines in Fig. 1). 

To study the spin torque on the FM moments, we pro- 
ceed in the following in two steps: First, we consider 
the TI without a ferromagnet, and briefly discuss the 
non-equilibrium spin density (5) ncq in the presence of a 
current in the TI surface (inverse spin-galvanic effect). 
In a second step, we use this spin density as a boundary 
condition and analyze the diffusion of the spin density 
into the ferromagnetic metal and the resulting torque on 
the FM moments. This two step approach is justified, as 
long as the chemical potential /irj 3> A cx , with A cx the 
exchange coupling of the surface-state spins with the FM 
moments [situation in Fig. 1(b)]. Note that the presence 
of spin diffusion makes the TI/FM bilayer different from 
a TI/FI structure. For the case of a TI/FI bilayer, the 
spin density merely acts as an effective field coupling to 
the magnetic moments of the ferromagnet which leads to 
a torque of the form T = A ex m x (5) neq , with m the mag- 
netization direction in the ferromagnet. This thus adds 
to the torque from the Oersted field (field-like torque) 
[6]. We will show in the following, that there is an ad- 
ditional perpendicular torque due to the spin diffusion 
(Slonczewski-like torque [18]). 

Dirac state - A TI surface state is described by the 
Dirac Hamiltonian 

^=W F (zX(T)-k-fl D (1) 

with a the Pauli matrices acting in spin space and z is the 
unit vector in z direction. This leads to the characteris- 
tic linear dispersion and spin helicity around the Fermi 
surface, see Fig. 1(c). The velocity operator v = d^'H^ 
is directly proportional to the spin operator 5 — (h/2)a 
and reads 

2 

v=-v F (zxS). (2) 
n 

While the TI has a vanishing equilibrium spin expec- 
tation, any current running through the system leads 
to a finite spin density independent of the microscopic 
details of the transport. We consider a current density 
j. x — en(v) ncq in x direction with e the electron's charge 
and n the electron density, which yields a spin density 

(Sy)neq = ~2eVF^' ^ 

If the magnetic moments of a FI couple to the surface- 
state spins through T-l ex = —A ex m ■ 5 as discussed in 
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FIG. 2. Spin accumulation in the ferromagnet (d = 8nm) as 
a function of distance z from the TI/FM boundary, where 
the solid (dashed) line denotes S± (S z ). For these plots, we 
used a spin decoherence length of A^ = lnm and the spin 
diffusion length of Permalloy A s f = 5nm[20]. Green (red) 
curves correspond to a spin-precession length Aj = lnm (Aj = 

0. 5nm), respectively. 

Yokoyama et al. [6], the spin polarization on the TI sur- 
face leads to a spin torque of the form T — A ex rhx (5) noq . 
This leads to a field-like torque, which for an in-plane 
magnetization is out-of-plane. In the next section, we 
show how spin diffusion into the FM additionally leads 
to an in-plane torque (Slonczewski-like torque), in a way 
similar to the spin-current injection in HM / FM bilay- 
ers [1-3]. 

Spin diffusion - Given the spin polarization at the TI 
surface, Eq. (3), as an input, we consider the diffusion of 
these spins into the ferromagnetic metal and the torque 
they thereby exert. The diffusion (in z direction) leads 
to a steady-state (itinerant) spin density through [19] 

= - V- Ji - — (5 x m)i - — [m x (5 x m)]* - — , (4) 

where the spin current (for the zth spin component) is 
given by 

Ji = -VVSi (5) 

with T> the diffusion coefficient. The second term in 
Eq. (4) describes the precession of the spins around the 
moments of the FM with tj the spin precession time. 
The third term captures the relaxation of the spin com- 
ponent perpendicular to m with the spin decoherence 
time, and the last term is the spin diffusion with time 
scale T s f. In the following, we use A s f = 5nm [20] and 
values for Aj and A^ of the order lnm (Af = X>t,). 

We solve equations (4) and (5) requiring no spin cur- 
rent through the outer boundary of the FM, J(d) = 0, 
where d is the thickness of the ferromagnetic layer. For 
the TI/FM interface, we assume that due to the exchange 
interaction, the itinerant spins of the FM right at the in- 
terface align with the spin density of the TI interface, 

1. e., 5(0) = 7(5) noq with 7 of order one [21]. With these 
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-Aj = lnm , A^ = lmn 
-Aj = 0.5nm , Xj, = lnm 




FIG. 3. (a) Total torque as a function of the FM metal thick- 
ness d. We used again A s f = 5nm, and the solid (dashed) 
lines denote the in-plane (out-of-plane) torque, (b) The ratio 
of the out-of-plane to the in-plane spin-torque efficiency as a 
function of X3/X4,. 



boundary conditions, the spin distribution in z direction 
is given by 



S(z) = S±(z) + iS z (z) = Sq 



cosh[fc(z — d)] 
cosh (fee?) 



with 



k = JX7 2 -iX 



J > 



(6) 



(7) 



and A^ 2 = X~ { 2 + X- 2 . S ± (z) is the in-plane spin density 
and So is the initial spin density [d = 0), both perpen- 
dicular to rfi. Figure 2 shows the in-plane spin density 
S± perpendicular to the magnetization (solid line) and 
S z along the z axis (dashed line) for d = 8nm. Note that 
this thickness d ~ 8nm 3> 1/k' with k = k' + ik". Using 
Eq. (6), we can thus approximate 

S(z) « S e~ kz = S Q cosk"ze- k ' z + iSosmk"ze- k ' z , (8) 

i.e., both components oscillate and decrease exponen- 
tially, see Figure 2. 

Figure 3(a) shows the torque as a function of the FM 
layer thickness d. The torque is given by the spatial 
change of the spin current compensated by the spin re- 
laxation, 



T = dz 
Jo 



- d z J(z) 



1 

r s f 



S(z) 



(9) 



Given the spin distribution in z direction of Eq. (6) , we 
find 

1 % , sinh(fccT) 



T = S ( 



A 2 
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k cosh(kd) 
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(10) 
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For the limit in the last line, we used d — > oo. As expected 
from the fast decay of the spin density in Figure 2, the 
torque is 'deposited' within only a few nanometers. The 
total torque exerted on the ferromagnet as a function of 
the thickness d thus stays constant with layer thickness. 

For the geometry described in Fig. 1(a), the spin po- 
larization perpendicular to the magnetization of the FM 
is y/2/2 of the total polarization (S y ) neq , and we find for 
the thick-FM limit (d > 1/k') 
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In analogy to the spin-Hall angle $sh = (2eJs)/(hJc), 
which describes the spin-Hall current per charge current, 
we define the spin-torque efficiency 
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jx h 
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(13) 



Figure 3(b) shows the ratio of out-of-plane to in-plane 
spin-torque efficiency as a function of Aj / A^ to illustrate 
that for Aj ~ A^, the two are roughly of the same size. 
Using Aj = A^ = lnm, A s f = 5nm, v-p = 5 x 10 ms _ , 
and a typical diffusion coefficient T> = 1 — 10cm 2 s , we 
find for the in-plane and out-of-plane-torque efficiency 
= 0.15 - 1.5 and \0 Z \ = 0.065 - 0.65. 
Rashba-split 2DEG- The real TI/FM interface not only 
hosts a topological Dirac state as considered above, but 
also a Rashba-split 2DEG with a dispersion 
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+ a{z x (?) ■ k — /ir. 



(14) 



Figure 1(b) shows schematically the combined bandstruc- 
ture with chemical potentials /xjj > //r > 0. Note that 
according to experiment, both vf and a are positive, 
such that an alternating spin structure with the spins 
clockwise on the TI surface state as denoted in Fig. 1(c) 
is obtained [11]. The spin-torque due to Rashba SOC in 
the presence of a current has been studied both theoreti- 
cally [22-24] and experimentally [25], and we use here the 
result for the non-equilibrium spin density in the strong- 
spin-orbit-coupling case [23], 



("SjJneq 



h maj x 



(15) 

2e 2E F v ; 

Together with the Dirac surface state, the total spin po- 
larization yields 
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Note that the contribution of the Dirac and the Rashba 
surface state are of opposite sign. Using a rough estimate 
for the current ratio for the two states j^/j^ ~ 0.6 [17], 
ma = 2(fc F - kp w 0.02A" 1 and E F w 0.2eV from 
ARPES[11], we find (2E F )(ma R ) w 5 x 10 5 m/s w u F . 
Both contributions are roughly of the same magnitude, 
and potentially cancel each other out. 

However, in as long as spin torque is measured, the 
dominant source between topological Dirac surface state 
and the Rashba-split 2DEG can be assessed through a 
direct measurement: The field- like torque coming from 
the spin polarization will either increase or decrease the 
field-like torque coming from the Oersted field. For 
Vf, ol > [11] and assuming a coupling through ferromag- 
netic exchange, the topological surface state adds to the 
torque generated by the Oersted field, while the Rashba 
surface states will weaken the torque. This thus allows 
for the identification of the TI surface state that domi- 
nates the generation of torque. 

Conclusions - In this work, we analyzed TI / FM bi- 
layers and their potential for magnetic-memory applica- 
tions based on large spin-torque generation. We pro- 
pose measuring the spin-torque efficiency of the TI sur- 
face states through ferromagnetic resonance as recently 
done for heavy metal/FM bilayers [1-3, 26]. Considering 
itinerant spins that diffuse into the FM, we predict both, 
an out-of-plane (field-like) torque, as well as an in-plane 
(Slonczewski-like) torque in such TI/FM bilayers. This 
is in contrast to what is expected of TI/FI bilayer (only 
field-like torque) and what is found in heavy metal/FM 
bilayers (primarily Slonczewski-like torque) [1-3, 26]. In 
our model, both components of the torque originate from 
the combination of the inverse spin-galvanic effect of the 
TI surface and spin diffusion into the FM. For metal- 
lic TIs, an additional in-plane torque could arise from 
the bulk spin-Hall effect. As transport is dominated by 
the surface states for thin TIs [17], we still expect the two 
components of the torque to be of comparable magnitude 
[see Figure 3(b)]. For realistic parameters, we predict a 
spin-torque efficiency of the order of \8\ « 0.1 — 1. This 
is comparable to or larger than the largest value of spin- 
torque efficiency observed to date [1-3, 26]. Therefore, 
our predictions point to an exciting prospect for TI/FM 
bilayers in magnetic memory applications. 

Finally, we comment on limits of the applicability of 
our approach to extremely thin FM layers and future 
directions. As the total spin torque stays constant in- 
dependent of FM layer thickness for d > 2nm, thin FM 
layers are preferable for device applications. However, 
our calculation treating the FM layer in z direction to be 
in the diffusive regime relies on a FM layer that is thicker 
than its diffusion length. For a device with an FM layer 
thinner than the diffusion length, the device should be 
modeled using a quantum Boltzmann approach [27, 28]. 
This will be a subject of future investigation. However, 
for the first experimental demonstration of the principle 



and the potential of the TI/FM bilayer, our simple model 
can guide ferromagnetic resonance measurements which 
do not require such thin FM layers. 
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